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Abstract 

Determination of all the coefficients in the crack tip field expansion for monoclinic 
materials under two-dimensional deformation is presented in this report. For monoclinic 
materials with a plane of material symmetry at *3 = 0 , the in-plane deformation is decoupled from 
the anti-plane deformation. In the case of in-plane deformation, utilizing conservation laws of 
elasticity and Betti’s reciprocal theorem, together with selected auxiliary fields, T-stress and 
third-order stress coefficients near the crack tip are evaluated first from path-independent line 
integrals. To determine the T-stress terms using the 7-integral and Betti’s reciprocal work 
theorem, auxiliary fields under a concentrated force and moment acting at the crack tip are used 
respectively. Through the use of the Stroh formalism in anisotropic elasticity, analytical 
expressions for all the coefficients including the stress intensity factors are derived in a compact 
form that has surprisingly simple structure in terms of one of the Bamett-Lothe tensors, L. The 
solution forms for degenerated materials, orthotropic, and isotropic materials are also presented. 

Introduction 

The use of fracture mechanics to assess the failure behavior in a flawed structure requires 
the identification of critical parameters which govern the severity of stress and deformation field 
in the vicinity of the flaw, and which can be evaluated using information obtained from the flaw 
geometry, loading, and material properties. In the linear elastic solids, stress intensity factors, k, (/ 
= I, n, HI), represent the leading singular terms in the Williams eigenfunction expansion series 
near a crack tip. k, are often assumed to be unique parameters associated with crack extension. 
The physical implications of the higher-order non-singular terms have been noted by Cotterell 
(1966). Especially, the so-called T-stress, second term of the crack tip stress field which 
represents the constant normal stress parallel to the crack surfaces, has been found as an 
additional parameter in characterizing the behavior of a crack (Larsson and Carlsson, 1973; Rice, 
1974). Cotterell and Rice (1980) showed that T-stress substantially influences the fracture path 
stability of a mode-I crack. The stress biaxiality parameter (Leevers and Radon, 1982; Sham, 
1991) has been tabulated as a function of relative crack lengths and overall geometry in many 
fracture test specimens for the isotropic solid using computational techniques (e.g., Kfouri, 1986 
and Sham, 1989 and 1991). Kardomateas et al. (1993) examined the third-term of the Williams 
solution and concluded its significance in the center-cracked and single-edge specimens with 
short crack lengths. 

In anisotropic linear elastic solids, Gao and Chiu (1992) examined the T-stress term of a 
crack in infinite orthotropic solids under mode-I loading. Because of the material anisotropy 
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involved, the T-stress term is affected by the material properties. It is also expected that, in 
general, mixed-mode crack behavior and the biaxiality parameter are also dependent on the 
material anisotropy. Thus, it is essential to develop efficient computational techniques to 
determine T-stress term coefficients including the stress intensity factors in anisotropic cracked 
materials with finite geometry. In this report, two methods based on the /-integral and Betti’s 
reciprocal theorem are proposed to obtain compact forms in calculating all the stress coefficient 
terms in the crack tip field expansion for monoclinic materials with a plane of material symmetry 
at *3 = 0. To determine T-stress term using the /-integral, the method by Kfouri (1986) is 
extended to anisotropic solids. The closed form solution of the auxiliary field, a point force 
acting at the crack tip, is derived for this purpose. A path-independent integral based on the 
Betti’s reciprocal work concept has been used for determining the stress intensity factors by 
Stem, Becker, and Dunham (1976), Hong and Stem (1978), Sinclair, Okajima, and Griffin 
(1984) for isotropic materials; Soni and Stem (1976) for orthotropic materials; and An (1987) for 
rectlinearly anisotropic materials. This path-independent line integral is also extended to 
determine all the stress coefficient terms with auxiliary fields. 

Mathematical Formulation 

In a fixed Cartesian coordinate system x,, (i = 1, 2, 3), consider a two-dimensional 
deformation of an anisotropic elastic body in which the deformation field is independent of the jc 3 
coordinate. In this report, attention focuses on the monoclinic material having three mutually 
perpendicular symmetry planes and one of the planes coinciding with the coordinate plane xj = 0. 
In this case, the in-plane and out-of-plane deformations are uncoupled. For in-plane deformation 
the strain and stress relations can be written as 

e=s'a (1) 

where e = [e, , £ 2 , y 12 f, <x = [<7„ , cr 22 , cr 12 ] 7 
or 


£, = s'j a j , i,j= 1,2,6 

where s' = s' are reduced compliance coefficients defined by s' = s y - s j3 s j3 / s 33 . 

Throughout the report, all indices range from 1 to 2 and the summation convention is 
applied to repeated Latin index unless otherwise noted. The bold-face letters are used to represent 
matrices or vectors. A comma stands for differentiation; overbar denotes complex conjugate. A 
symbol Re stands for real part; Im for imaginary part. 

In the absence of body forces, general solutions of the displacement vector u, the stress 
function 0, and stresses a, for in-plane deformation, according to Stroh formalism (Ting, 1996), 
can be represented by 
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( 2 ) 


u = Re[^a a d a f(z a )] 
0=1 

0 = Stliwwi 

a=I 


or 

u = ReU(/(z)}d] 

« = Re[B(/(z)>rf] 

Ojj — ~(p, y 2 ’ ^i2 ~ $',l (^) 

where 

(/(Z>) = rf»flg[/(z,), /(z 2 )] 

= *, + P a x 2 ’ JmtMa l > 0 

f(z) is an arbitrary function, d is a unknown complex constant vector to be determined. fl a , 
a a , and b a are the Stroh eigenvalues and corresponding eigenvectors determined by elastic 
constants only. For in-plane deformation, [l a are given by the roots of the characteristics 
equation: 

4a* 4 - 24A* 3 + (24 + 4,y - 24 a* + 4 = 0 (5) 


with positive imaginary parts. From energy consideration, Lekhnitskii (1963) showed that the 
roots are either complex or purely imaginary and cannot be real. A and B are Stroh matrices 
given by 




P\ 

<7i 


Pi 

Pi 


-ft 

-ft" 

B ] - 1 

-1 -ft 

1 

1 

> c — 

ft -Pi 

. 1 Pi . 


Pa = 4 lA*„ 2 - 4A*c + 4 . Pa = 4ft* ~ 4 + 4 1 Pa • 


( 6 ) 

( 7 ) 

(S) 


The eigenvectors a a and b a are unique to an arbitrary multiplier. Introducing normalization 
factors k a , we have 


A = 

^\P\ ^2 P 2 

, B = 

KP\ 

1 

K> 

3= 

to 


^ 2 



k 2 j 


The values of k a satisfy the conditions 

bp = 5 a p 
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Based on the Stroh formalism, the matrices A and B defined in eq. (9) satisfy orthogonality 
relations. For in-plane deformation, these relations can be expressed by 

B t A + A t B = I (10) 

B t A+A t B= 0 (11) 

where 


It can be proved that 


2 &,"(<?, 1, 2k;(q 2 1 

B~ T B~' = -2 /IT 1 , L ' =-Im[AB'] 


( 12 ) 


where L is a real, symmetric, and positive definite matrix which will be used frequently in the 
sequel. 

Note that the normalization factors cancel each other for the term AB l in the second 
equation of (12). Therefore there is no need to introduce the normalization factors in computing 
L 1 . From eq. (12) 2 , it is easy to get 


L' 




1 


s' ( be-d 2 ) 


e 

-d 


-d 

b 


(13) 


jl^+ H 2 = a + ib, fl { j± 2 = c + i d , 


e- ad- be- Im + i u i)] 


(14) 


For crack problems, it may be convenient to introduce a complex potential function <£ 
(Guo, 1991) such that 


<b = B(f(z))B-'g (15) 

where g = Bd, then the displacement expression in eq. (3) and stresses in eq. (4) can be rewritten 
in terms of the potential function O 

u = Re[Afl-'0] (16) 

=- Re[ < J > ,, 2 ]. cr, 2 = Re[0 (1 ] (17) 


The traction vector t at a point on a curve T with unit outward normal n is given by 


t . = - Re 


, t = - Re 

~d<f>~ 

i 

ds 


ds 


(18) 


where s is an arc length measured along F as shown in Fig. 1. Thus, without loss of generality, 
the traction free boundary conditions on a boundary may be written as 


RefO] = 0 on r 


(19) 




Fig. 1 The surface traction tona curved boundary r with a unit outward normal vector n. 

The resultant force and moment about the xj axis due to the surface traction t acting on T 
between Si and s 2 (s 2 > si) are 


[ ' t(s)ds = Re[<J>(s.) -<X>(5,)] 

(20) 

J '(*/-> - x 2 t x )ds = -Re[x,<I> 2 ~xT' 

(21) 

where 


X(Z) = fo 2 a)dX 


If r encloses a region and there are concentrated force/ and moment M inside the region, then 
the equilibrium of the body demands that 

ii 

% 

, u 

1 

(22) 

- ( jv, - x 2 t x )ds = M 

(23) 

Crack-tip fields 


Consider a crack in the anisotropic body. Let a coordinate system be attached to the crack 
tip and crack plane lies on the xj - X3 coordinate plane. The configuration is shown in Fig. 2. The 
crack faces are assumed to be traction-free. Note that the crack plane may not coincide with the 


symmetry plane of the material. 
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X- 

crack 


Fig. 2 A cracked body and a contour around a crack 


To find the solution for the crack tip field, employing eq. (15) and (16) with 
f a (z a ) = Za +I , we take the solutions u and O in the form 

« = Re[A(z 5+l )£?-'£] (24) 

O = b{z s+ ')B' 1 g (25) 

where the complex variable z a is defined by 

z a = r(cos6 + /i a sin#) , -k<Q<k 

g is a complex constant vector and <5 is the complex constant. We seek admissible values of <5 
subjected to a restriction in which the strain energy is bounded as r — » 0, that is 

- 1 < Re($) 

To maintain a unique solution for <I> in eq. (25), we introduce a branch cut along the 
negative x, axis with a convention, that is 

z s a + ' = r s "e ±i{S+m at 9 = ±k 

Insertion of eq. (25) into traction-free boundary conditions 

Re[$] = 0 at 6 - ±n 


yields 


r S-5 e ^ S+ 5 + 2 )g+ - = 0 
r S-5 e -^5 + 5 + 2) g + g =Q 


Subtracting of eq. (27) from eq. (26) leads to 

[^i2(<5+5+2)ff _ J]Jg = o 

where I is a 2 x 2 identity matrix. For a nontrivial solution of g, we must have 

det[(e i2(5+f+2) ' r - 1)1] = L' 2 <^ +2) * . i)i| = o 
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or 


Re(5) = (n-2) / 2 , n = 1,2,3,- 

(30) 

With eq. (30), from eq. (26) or (27), we have 


8 

_ _ r , 2 inx<5) C0S [27rRe(<5)]# 

(31) 

Since # is a constant, from eq. (31), (30), 


Im(<5) = 0 , 8 = (n-2)/ 2 

(32) 

g = < 

! o 113 

8 ' S ~ 2 ’ 2’ 2 ’ 

[-#, 5=0,1, 2, •• 



that is, 8 is real, the vector# associated with 8 is real if <5= - 1 / 2 , Vi, 3 / 2 , ...; and# is pure 
imaginary if <5 = 0, 1, 2, .... It is clear from the determinant given in eq. (29) that each of the 
eigenvalues 5 is a root of multiplicity two. Since # has two arbitrary components, we have two 
independent eigenfunctions associated with the double eigenvalues 8. Therefore the assumed 
forms of u and O given by eq. (24) and (25) are justified; no logarithmic type of solution form 
exists. 


Superimposing all the solutions with different orders of r, the crack tip field can be 
constructed as 


m = ^ R e[A(z 5 " +l )j? _, #J 

n= 1 

4> = '£b{z‘--')b-' s , 


(33) 


where 8 n = (n-2)/2.g n is real for n = 1,3, 5, — , #„ is pure imaginary for n = 2, 4, 6,—, and #„ 
are dependent on the geometry of the cracked body, material properties, and loading conditions. 
Similarly, there is no need to introduce k a in calculating Cfy and m, from eq. (33). 

Performing algebraic calculation and defining g n = [# nl , # n2 ] T , the stress and 
displacement components can be written as 

<7.! =X& +1 )' A Re { IsnMd" 20.02 farf' -A^f")]} 


n= 1 




22 = X ( s « + 0 Re { Isn 1 ki - S?" )+ 8 m farf* - fhtf )]} 

n=I (.A 1*2 J 


(34) 


C,2 = X ( S n + 0 r S " Re i L )+ *,2^2 fef" - d" )]| 




.A*. -1*2 


1 



u, = V r' 5,,+I Re< 

{Pi~Pz 


n=\ 


« 2 = X r<5,+l 


1 


“ {Pi ~P 2 


■ t,l (P2?2 n+1 " P,?f” +I )+ gnlfPxPz^ ~ PzPl^ )]| 

b.l fc?2” + ‘ - <7lSf” +1 )+ - P 2 ^f" + ' )]| 


where g a = cos 6 + pi a sin G . 


( 35 ) 


For the leading-order term, dj = -1/2, letting g, = -> JlUt k , the singular solution and T- 
stress term for the monoclinic solids can be written in the form: 


<7 = 






+ 0(r U2 ) 


(36) 


+ k„£ { ,' l ) )+T 

•sJ27tr 


3 12 


M6 


+ 0(r vl ) 


(37) 




p<» = s 'a m 

c n — ■> U // 


W = 


= J—fc/®/” +*//»//' ) +7> 

K 


cos 0 + 5 ,' 6 sin 6 
sin 6 


+ Q)r 


-sin0 

COS0 


+ 0(r vl ) 


(38) 


where cr*"’ , « (n) ,■ ■ -,are the /z-th order terms of the crack tip field; the subscripts I or II indicates 
the distribution associated with mode-I or mode-H k = [k lf k, ] r ; ki and kn are stress intensity 
factors for mode-I and mode-II respectively, ft) is a constant representing rigid body rotation, and 


(<T„), =Re 

PiPz 

( 

Pi 

V 

Pi 

(^n)//= Re 

1 

( jL. 

2 M 
Pi 

Pi -p 2 



Px-Pz 

\~J$2 


(<7 22 ), = Re 

1 

( 

Pi 

Pz ' 

, (^ 22 )// = 

1 



Pi -p 2 


"^1 

Pi -p 2 



§ 

to 

II 

50 

rt» 

P 1 P 2 1 

1 

__l_ v 

(^ 12 )// = 

1 

f 

JP J 

v 

Pz 

1 

J: 

l 

“*5 

w 


V£T,. 

Pl-Pz 


Jtz [ 
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(«,), = Re 

_ Ml M2 

(^i ) // 


(w 2 ), = Re 

iV^T) 

LjU : -/X 2 'J 

(^2 )// 

ifi 2 Qx'joi) 

J 


For the second-order term, & = 0, 


u 


( 2 ) 



cos 0 + s' 6 sin 6 
■y,' 2 sin 9 


+ (dr 


-sin# 

cos# 


(39) 


0-, ( , 2 } = T , <J { £ = cr, ( 2 2) = 0 (40) 

where a represents the rigid body rotation, and 

T = -i(g 2 ib + g 22 d), (0 = -is' n (g 2i d + g 22 e) (41) 

Clearly the T-stress term is dependent on the material properties in anisotropic solids. 

For orthotropic materials where xi and X 2 coincide with the material symmetry axes, the T- 
terms become 


w (2) = Tr 


s' l{ cos# 
s,' 2 sin 0 


(Or 


-sin# 

cos# 


a, ( f ) = r 


In isotropic solids, the second-order terms are 


( 2 ) 

u = 


Tr 


cos# 
-v sin# 


+ (0r 


-sin# 

cos# 


, (T| 2) = T for plane stress 


u a> _ nuvy 


(l - v)cos# 
-vsin# 


+ 0)r 


-sin# 

cos# 


, 0'| < 1 2) = T for plane strain 


For the degenerate materials, the crack tip fields are derived in Appendix A. 


Auxiliary Fields associated with the Crack Tip Fields 

Some auxiliary fields with higher-order singularities are needed in order to determine the 
coefficients in the expansion of the crack tip field by the use of conservation laws of elasticity 
and Betti’s reciprocal theorem. Since the negative integers are also the eigenvalues of the crack 
problem which satisfy zero traction on the crack surfaces and satisfy the field governing 
equations of the anisotropic solids, the associated eigenfunctions can be conveniently used as 
auxiliary (pseduo) fields. Note that each eigenfunction with higher-order singularity has 
unbounded energy near the crack tip and thus corresponds to some concentrated source at the tip. 
This eigenfunction can be imagined as a self-equilibrated solution to the crack problem under 
some specified loads. These auxiliary fields may be obtained by choosing the values of n in eq. 
(33) as negative integers, that is, 


9 



u‘ =Re[A(z 1 -* l )B‘ , *»] 


\ / 

=B(z Am * x )B~'h m 

(42) 

= -Re[0>, 2 ], < 7 “ = Re[0, ,] 

(43) 

A„, = ~m! 2,m= 1,3,4,-, 

(44) 


K=IK,,KJ 


where h m are arbitrary constant vectors. h m is real for m = 1, 3, 5, •••; h m is pure imaginary for m = 
4, 6, Utilizing eq. (22-23), the auxiliary fields, except m = 4, defined by eq. (42) yield zero 
resultant force on any contour 7~ which encloses the crack tip shown in Fig. 2. The corresponding 
resultant moment about the Xj-axis, produced by the tractions acting on the contour r is also zero 
for the auxiliary fields in eq. (42). The special case of A n = -2 or m - 4 can be also directly 
explained from eq. (23). In this case, the function associated with [0, /i 4 , ] r corresponds to the 
particular solution for a crack under a concentrated moment about x 3 -axis, (-2 n i h 42 ) , applied at 
the crack tip; the function associated with [h 4l , 0] r corresponds to the homogeneous solution 
which satisfies zero concentrated force and moment at the crack tip. 


From eq. (42)-(44), the stress and displacement components of the auxiliary fields are 


= (A, + l)r A " ReJ — - — [h mX (p i;g A " - rfgf * )+ h m2 ^ 2 (/i 2 g A ” - )]1 


. Ah 1^2 


=(A m + l)r A " Re< - ~ [h mX (g^ -g^ )+ h m2 (pi x g^ ~n 2 g x " )]| 

=(A m + l)r A » R e|— ^ — [h mX (jx x g^ - n 2 g *"' )+ h m2 n x p 2 (g*“ - gf* )]J 

u‘i = r A ’"+' Rej— 3— k.U^ A "' +1 -P^ m+ ')+h m2 (M x P 2 <;2 m+] 

Mi ^2 1 


(45) 


a* = r ^ +1 R e J 


1 


Ml-»2 


kl k?2 m+I - )+ Kl (^.^2 A " +1 - )][ (46) 


In the following two sections, stress intensity factors, T-stress term, and coefficients of 
higher-order terms are determined using the 7-integral and the Betti’s reciprocal theorem 
including the use of above auxiliary fields. 


J-Integral 


(a) T-stress Term 



/ conservation laws (Knowles and Sternberg, 1972) for a plane anisotropic elasticity 
problem may be written as 

J k = J (wn k - t,u i<k )ds = 0, it— 1, 2 (47) 

for an arbitrary closed contour C that encloses no defects, cracks, or material inhomogeneities. In 
the above equations, W is the strain energy density, W = cr, ; £, ; / 2 , where and £ :J are the stresses 

and strains respectively; t, are the traction components defined along the contour, t, = <7,, n t ; n k are 
the unit outward vector normal to the contour path. Letting k = J, the conservation law is reduced 
to the Rice’s path-independent /-integral or the rate of energy release rate per unit of crack 
extension along the xi - axis, which is given by 

/ = J r (<7 r £n, 12 - 1 T u (48) 

where F is an arbitrary path which starts on the straight lower face of the crack, encloses the 
crack tip and ends on the upper straight face with the positive direction in a counterclockwise 
direction shown in Fig. 2. Here, the crack surfaces are assumed to be traction-free. 

Consider a cracked body under the two-dimensional deformation. The components of 
stress, strain, and displacement are represented by <7 iJt % and u„ respectively. As r — > 0, the 
asymptotic fields including the constant T-stress terms are given before. Now the coefficients of 
the T-stress terms and third-order terms are derived using conservation laws. 

In general, for the purpose of determining the coefficients g„ of the term r s " (8 n > -1 / 2) 
in the actual crack tip stress field, we may employ the /-integral method and follow the following 
procedure: 

(i) find an auxiliary (pseudo) field that has singularity ~ r -5 " -1 as r — >0. It is convenient to 
select auxiliary stress field which gives zero traction on the crack surfaces and contains only 
the stress singular term r ~ s "~ l ; 

(ii) superimpose the actual field (the mixed-mode boundary value problem, in general) on the 
auxiliary field and represent the /-integral for the superimposed state as 

J, = ]+J, + J„ (49) 

where 

J s = [(cr +<T a ) r (£ +£ fl )«, /2-(t + t a ) T (u i +«“)]/5 

J a =\ r [{o a ) T e a nJ 2-(t a ) T u^]ds 
and 

/* =/,-/-/. 

= {[(T r £°+((T fl ) r £]n 1 /2-t T U°,-(t a ) T U ] }ds 
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(50) 


= J r [( tVw, -t T u‘ -(t a ) T u A ]ds 

= \ r ( <y l jKj n \- t K\- t “ u i.\) d s 

where the superscript or subscript “a” denote quantities referred to the auxiliary field; / is the J- 
integral for the superimposed state; J for the actual state; and J a for the auxiliary field and J M is 
the interaction integral. In the sequel, we assume that the J-integral is path-independent for both 
the actual field and the selected auxiliary fields, denoted by J and J a . Then the integral J s for the 
superimposed state, thus Jm, is also path-independent. If the auxiliary fields given by eq. (42) are 
used, it is readily proved that 

7*0, for A„ =-1/2 

f.-O. for A < 5I > 

(iii) evaluate J M as T — >0. For simplicity, r may be taken as a circle with radius r, as r — > 0, the 
only terms in the integrand that contribute to J M are the cross terms between r 5 " in the actual 
stress field and the auxiliary stress term with order r~ 5 " _l ; 

(iv) carry out the routine manipulation, the exact expression for J M can be obtained as 
Jm =J M (gn ) when r — » 0; 

(v) evaluate J M for a finite contour r using the computed actual field and the exact auxiliary 
solution; and determine the coefficients g n from the value of J M and the expression of 

J M = J M (gn ) SS r — > 0. 


In extracting the T-stress in eq. (33) org 2 , we make use of another auxiliary field, that is the 
solution to a point force /(per unit thickness) applied at the crack tip. Under the load, cr° « r' 1 . 

Note that the point force / must be resisted by traction t applied to some boundary C in achieving 
equilibrium. In the Stroh formalism (Ting, 1996), the real form solutions due to the point-force 
application can be written as 


2 u a 


Inr 

K 


1 + 5(0) 


(52) 


20° =L(9)h 


(53) 


where h = L' 1 / , / = [/,, / 2 f 

S(6) = — Re[A{ln(cos0 + /isin0)Z? r ] 
L{6) = - J- Re[£{ln(cos d + sin 6)B t ] 


X] = rcosd, X 2 = rsind. 
It assumes the values 
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0 0=0 

ln(cos0 + p a sin0) =1 . 

± in 6 =±n 


In a cylindrical coordinate system (r, 0, xy), let t r and t$ be the traction vectors on a 
cylindrical surface r = constant and on a radial plane 0 = constant, then 

t r =-~K , t e =f r (54) 
r 

e° =n T t r , o a re =m T t r =n T t g , cr a e =m T t g 
where n T = [cos0,sin0] and m T = [-sin0, cos0]. 

It follows from eq. (52) and (53) that 

h=tr= 0 - = a °re =° 

(T® = ~n T <j>^ / r 


or 


/, Im f [(/ti^ 2 -l)sin0 + Oi, + J u 2 )cos0]cos 2 0-sin0 
iTir [ g,g 2 


f [(/r, n-, -l)cos0-(jU, + jU,)sin0lsin 2 0 + //,/x 2 cos0 
— — w : ► 

~n r SiSi 


(55) 


and 


u = — 



5 n In r 


5,, In r 


b 

d 

d~ 

e 


- Im 


- Im 


Pi I" Si ~ Pi ln Si 

[Pi - Pi L $2 ln S 2 - <?. ln Si JJ 

Pi Pi In g 2 - PiP I ln Si 
Pi~ Pi I PiSi In - P 2 q, In g x 


For isotropic materials, 

o a r8 = (7g = °, o “= — — (/i cos0 + f 2 sin0) 

nr 


tt a = 1 < 

8 7i G 

(K" + l)lnr 

-(*-1)0' 


-2 sin 2 0 

sin 20 

17 ," 

. Of -1)0 

(/c + l)lnr 


sin 20 

2 sin 2 0 J 

hi. 


(56) 


(57) 

(58) 


K = 3 - 4v 
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By superimposing the actual field on the auxiliary field, and using the path-independent 7- 
integral for the monoclinic elastic cracked body, it can be proved that 

J s = J + T s' u /, + cof 2 

J M =Ts' u f ]+ C0f 2 (59) 

Note that because a i 0 ' r " 1 

K =J r -»X.)*=° 


Kfouri (1986) used the method to calculate the T-term for isotropic materials. Wang et al. 
(1980) and Wu (1989) applied the 7-integral to determine the stress intensity factors for 
rectilinear anisotropic solids and general anisotropic materials respectively. In this report, the 
method is extended to determine all the coefficients in the crack tip field expansion for 
monoclinic materials. From eq. (59), it follows that 


T = 



S'u /, 


and 


co = 


u 


Detailed proof of eq. (59) is given below: 

For general anisotropic linear elastic solids, we have the following relations 


O £° = c u£„S a = <r“e and oe a =Ou a 

ij ij ^ ijkl kl ij ij ij ij ij tj t’J 


where c ijk , = c UlJ . From eq. (50), 


du° 


J u = l r < J «i - t, <i - «,, 1 ) < ds = J r (<7,2 - 1 ■ M, j ) ds 


(60) 


where du“ / ds is the tangential derivative of u 1 ’ . As r —> 0, we evaluate J M and note that the 
only terms that contribute to J M are the cross terms between T and f. After substituting these 
fields into the integral and performing the routine algebra, the integral J M may be evaluated as 


r uU 

^ M = r^oJr (<7 ' 2 "^9 t"u iX )ds 

=J r « 


ds 

^^L_ t a u Cl) )ds 


ds 

=-J r ds 

= -J r («,| 2) ) r ^' ds 
= -( M f) r J r r ds 


= (U?) T f 


( 61 ) 


In the above derivation, eq. (22) has been used. From eq. (33) and (39) via ( 12 ) 2 , 
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(62) 


« 


( 2 ) 


= -iL~' g 2 = 


Ts' 


0) 


Insertion eq. (62) 

into (61) leads to 



J M =-ig T 2 L-'f = -if T L' ] g 2 

(63) 

and 

J M = Ts' u f l +Q)f 2 

(64) 


Eq. (63) and (64) can be used to calculate g 2 and T. As/ is arbitrary, it is convenient to 
choose/ to be the following values 

[1 , 0] T = e x , [0,\] T =e 2 

respectively. Note that e * has a dimension force/length. Correspondingly, eq. (63) yields two 
linear equations and they are, in matrix notation, 

K = L g 2 

where 

7 = r/ (,) y (2) i r 

J M L J M > J M J 

and J is the value of J M when / = e*. Therefore, 

g 2 = iLJ M (65) 

Using eq. (64), The T-stress and co can be obtained as 

T = J™/s' u , c o = J ™ ( 66 ) 

However, the choice of auxiliary field is not unique. For instance, we may choose the 
auxiliary fields as the sum of the field (<7 ',m') for the point force/ and a field (cr ",u") which is 
a known solution for the same cracked body under some loads on the outer boundary, then 

o a =<7' + <j" 
u a = u' + u" 

Superimposing the actual field on the new auxiliary field, we have 

J,[o + o a ] = J[o] +J a [cn+J M [c,c a ] (67) 

From the definition of J M , 

J M [<7,<j a ] = J M [(7,0' + a"] = J M [a,a']+J M [o, a") (68) 

J M [o ,c"]= J s [(7 + (7"] - J[<7] - J[<7 "] (69) 

In elastic materials, J is equal to the energy release rate G in the absence of body forces and 
dislocations and is related to the stress intensity factors (see Appendices B and C) through 


15 



( 70 ) 


G = —k T L~ ] k 
2 

Therefore, eq. (69) leads to 

J M [o,a"] = k T L~ ] k" (71) 

where k is the stress intensity factor for the state ( a,u ) , and k" for (cr", u") . 

Using eq. (71) and (64), eq. (68) can be expressed by 

J M [0,0 a ] = T s[J x +0)f 2 +k T L~ l k' (72) 

Inserting eq. (72) into (67) and utilizing eq. (13) yield 

J s [o + (j a ] = 7[cr] + / n [<7 n ] + T 5,', /, +cof 2 

+ Im[(//, + [i 2 )k 2/ k a + fi^fi 2 i,k uk j + + fl^k/k^ 

For isotropic case, 

J s [<y + e ° ] = J[o] +J a [o a ] + T s' u /, + O) f 2 + +2 s' u (k„k" + k,k',') (74) 

For plane strain under mode-I loading, if/= (/), 0), then eq. (64) and (74) reduce to the results 
given by Kfouri (1986) and j', =(1 -v l )!E 

(b) The third term 

The third-term coefficient in the asymptotic solution can be also obtained from the J- 
integral method. An auxiliary field with singularity <j“ ~ 0{r~ V2 ) can be introduced by selecting 

m = 3 in eq. (42). By superimposing the actual field (the mixed mode boundary value problem) 
on the auxiliary field, the interaction integral J M may be evaluated as 

J M =-^hlL~ ] g 3 (75) 

Following in a similar manner, 

g, =~Lj u (76) 

3 K 

where J M = [J™ , J (2) ] T and J™ is the value of Jm when h 3 =e k ,(k = 1,2)- e k has a dimension 
of force/(length) 1/2 . 

In general, superimposing of an auxiliary field with r A " on the actual field and 

applying the 7-integral to this combined state with derivations proved in Appendix D for n * 2, 
we can get the interaction integral J M denoted by J M „, that is 

J Mn =-2x8 n (d n +l)h T a L- } g n , n= 1,3, 4,5,-. (77) 
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Then 


where 


gn 


LJ Mn 

~2x8,J.3, + iy 
iLJ Mn 

„ In 5 n (S n + 1) ’ 


n = 1, 3,5, ••• 
n = 4,6,8, - 


/ =ry (1) 7 (2) i r 

J Mn \- J Mn ' J Mn J 


and J^l i s the value of J M „ when 


(78) 


f e k , n = 1,3,5,- •• 

| ie k , n = 4,6,8,--- 

Here, ey possesses dimension force / (length) 1 " 5 " . For the first singular term, introducing stress 
intensity factors, k = [k„,k, ] r =4nT2g x for the actual field andfc" = [k°, ,k°] T = 4x12 h, for 
the auxiliary field, J MI and k can be rewritten from eq. (77), (78) as 

J m =(k a ) T L' l k 


k = Lj m 

where j M] = [Jy\ , J ( 4\] T an ^ n is the value of J m when k a =e k . 

Betti’s Reciprocal Theorem 


(79) 


(a) T-stress Term 

For a linear elastic plane problem, Betti’s reciprocal theorem can be stated as 

J c (f ■ u a — 1“ ■ u)ds = 0 (80) 

where C is an any closed contour enclosing a simple connected region in the elastic body; u is the 
displacement vector and t the traction on C corresponding to the solution of any particular elastic 
boundary value problem for the elastic body; u a and f are corresponding quantities of the 
solution of any other problem for the body. Considering a crack in an anisotropic linear elastic 
material, and suppose the crack surfaces are free of tractions for both elastic states. If the closed 
contour C encloses the crack tip and extends along the crack surfaces, then it can be shown that 
the integral 

/ = j r (f • u a - t a ■ u)ds (81) 

is path independent where F is an any path which starts from the lower crack face and ends on 
the upper. Let (t, u) be an actual state for the crack under consideration, then eq. (81) provides 
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sufficient information for determining the amplitude for each term in the asymptotic crack-tip 
fields if proper auxiliary solutions (f, « a ) are provided. In this section the Betti’s reciprocal work 
contour integral is used for computing stress intensity factors, T-stress and other higher-order 
coefficients for monoclinic materials. The procedure can be evaluated from the analysis as 
follows. 

For determining the coefficients g„of the term r 5 n (£ n > -1 / 2) in the actual crack tip stress 
field, an auxiliary (pseudo) field with o° r ~ s ' : ' 2 or w" «= r~ 5n ~' can be chosen. As r — >0, take a 

F as a circle around the crack tip and evaluate integral 7. When r — >0, the only product between 
g„ and the auxiliary terms in the integrand given above can contribute to the integral 7. Therefore, 
the expression for 7 = I(g„) can be obtained as r — » 0. The value of 7 for a finite contour F shown 
in Fig. 2 is available from the numerical solutions for t and u of the boundary value problems and 
the exact auxiliary solution. The g„ can be computed from the expression for I = I(g„) and the 
value of 7. 

To determine the T-stress org 2 for the crack-tip field from eq. (34), the auxiliary elastic 
field with stress singularity cr" oc r~ 2 as r — > 0 is used and can be obtained from the eq. (42) by 
choosing m = 4, that is, in Stroh formalism, 

u a = Re[A( z - ] )B-\] 

/ , < 82 > 
=B(z~ ] }B- l h 4 

The moment about x,-axis applied at the crack tip, using eq. (23) and (82), is given by 

M = -2k i h 42 

When F shrinks to the crack tip, it is clear that only those parts of the integrand in eq. (81) which 
behaves like O (1/r) as r — » 0 can contribute this portion of the integral. Substituting these fields 
of the two states into eq. (81), performing the integration for the circle surrounding the crack tip 
and evaluating the results in the limit of vanishing radius, the results may be derived, and 

/ = lim f (f u a -t a ■ u\ds = -2k h T 4 V l g, (83) 

r— »o Jr 




(84) 


where 7 = [7 (1) , / (2) ] and I {k) is the value of I when h 4 =ie k ,(k- 1, 2). (Dimension of ey is 
force). 


From eq. (13), (41), and (84), 


7 (1) 
2 k s, | 


/( 2 ) 

CO 

2k 


(85) 


For isotropic materials, the auxiliary displacement vector and stress functions can be modified as 



<D 


( 86 ) 


sin0 1 + cos 20 
r sin 20 


sin 20 
1 - cos 20 


Im[/t 4 ] 


— a 

1 

2 

- 2 cos 30 cos 0 

-sin 40 

-.a 
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-2 sin 30 sin 0 2(cos20 - 1) sin 20 

a 

r 

-sin 40 

- 2 sin 36 sin 6 




(87) 



i 

- (k + 1) cos 0 + 2 sin 20 sin 0 

- (k: - 1) sin # - cos 20 sin 0 

u\ 

4 Gr 

( K - 1) sin 0 - cos 20 sin 0 

— (/C + l)cos0 -2 sin 20 sin# 


Im [/? 4 ] 


( 88 ) 


Then the path-independent integral I has the same form as eq. (83). Eqs. (84) and (85) are still 
valid. 


(b) The Third Term 

The coefficients of the third term in the eigenfunction expansion of the stress field can also 
be obtained using Betti’s theorem. Selecting m = 5 in eq. (42), an auxiliary field with stress 
singularity cr" ~ r~ 5 ' 2 desired for this purpose can be obtained. Applying the Betti theorem of 

reciprocity to the actual field and the auxiliary field and evaluating the integral / as T — » 0 near 
the crack tip, we obtain 

I = -3k h T 5 U'g 2 (89) 

Eqs. (89) will be used to calculate g 3 for mixed-mode problem when the two proper auxiliary 
field solutions are provided, g 3 can be expressed in the form 

= (90) 

3k 

Applying Betti’s reciprocal theorem to the actual fields and auxiliary fields with 
<t" « r A "* 2 , the path independent 7 denoted by can be evaluated by 


1 . .2 - - 2 * + 1 ) hi, £■' g. 


It follows from eq. (91) 


Sn = 


LI 


<1+2 


27r(<5 n +l) 
ijJj* 2 
2k(5„ + 1)’ 


where 


n — 1 , 3 , 5 , 
n = 2,4,6, 


(2)iT 


/,=[/, ( \ iy'] 


(91) 


(92) 


and Ij k) is the value of 7 y when 
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j e k , n= 1,3,5,--- 
\ie k , n = 2,4, 6,--- 


Here, possesses dimension force x (length) 5 " . The detailed proof is shown in Appendix E. 
For the first singular term from eq. (91) and (92), I 3 and k can be written as 

/ 3 = -7T hi L 'g , = -J2n hi L' 1 A: (93) 


k = - 




(94) 
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Appendix A 

Deformation Field for Degenerate Materials 

The solutions for the nongenerate materials can be modified so that they can be applied for 
degenerate materials. We write the general solutions for nongenerate materials as 

u = RefAZT'Xfzjg] 

O = X(z)g 

where 


»-i 


X(z) = B{f(z))B- 
In the limit /t, = jj^ = ji , it can be proved that the matrix reduces to 

X(z) = f(z)I + x 2 f\z)V 

where z = x t + jix 2 \ 


V = 


\i ll~ 
-1 -// 


Hence, for isotropic materials, we can obtain 

1 


u = ~ Re{[/(z)£ - ix 2 f'(z)v]g} 
0=[f(z)l + x 2 f'(z)v]g 


where 


E = 


S 11 ^12 


i2s n 


^11 ^12 


(*^1 1 S\2 ) l2.S\ j 


V = 


l 

-1 


-1 
— i 


For a crack in isotropic materials, choosing/fzj = z S+l and performing routine manipulations, the 
crack tip fields can be represented as 


'22 


<7,2 


= X(<5„+D^Re 


2e ns n e-K/2) + sin Q e ,6,fi + sin Q 

-8e i(S ”- m sin# e lS ’ e -de‘ U,2+(S ^ )e] sin# 


,[ff/2+(5„-i)e] , 


iS.e 


+ Se‘ 


[jr/2+(<5„-i)e] 


sin# 


■8e' 


(S n - D8 


sin# 


§ n I 

§, i2, 


r s.*' 

= Y — Re< 
“ 4 G 


' (2C + H_ 2(5„ sin0 (K-l)e ,iS -*' )e + 2{8 n + s in0 

+2{8 n + sin0 (ic + l)<r ,,r,:M *"* w -2(5, +\)e iS ’ e sin 6 


Sn I 

L«.:J 


where k = 3 - 4v for plane strain; K = (3 - v) / (1 + v) for plane stress. 
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Appendix B 
J-integral 

The /-integral with a contour T around a crack shown in Fig. 2 can be written as 


J=j r (W n,-f, «,.,)<& 



If/ du, . , 

= 2lr^l7" lUl ' )JS 

= ~j r (*I du-t T it , ) ds 

(Bl) 

Substituting the expressions fort, = 0 , u and t and using the orthogonality relations 


B t A + A t B = I, B T A + A t B =0 

(B2) 

J can be further expressed by 


J = — Re 
4 

+ •)(«. + I)//-* 8 '*}*'' 

_ m n 

(B3) 

It can be readily shown that 


l z 

*•**•*„ -ln(ij/ii) + i 2 * «.=«.= — 

(B4) 

The term contributing to J is the first term and 


/ = 77 Re{[ln(r, !r { ) + i 2n]gJ B T B ] gl } 
16 

(B5) 

By the use of the identity 



B T B l =- 2/ IT 1 

(B 6 ) 

eq. (B5) becomes 



J = 7 Im{[ln(r 2 / r, ) + / 2 n]g T x IT 1 g[ }= -g] ZT 1 g] 
o 4 

(B7) 

The above derivation has used the real values ofgi and L' 1 . Since 



Isi) ^ 

II 

bo 

(B 8 ) 

we have 

/=-* r r'A: 

2 

(B9) 
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Appendix C 
Energy Release Rate G 



Fig. 3 Crack opening displacements for the extended crack and stress distribution 
ahead of the crack tip prior to crack extension. 

Irwin (1957) realized that if a crack extends by a small amount A a , the energy absorbed in 
the process is equal to the work required to close the crack to its original length. Using a polar 
coordinate system with the origin at the extended crack tip, the energy released for a unit area to 
extend is written as 

1 fifl 

G = lim G i2 (Aa- r,0)[u j (r,7t)-uXr,-7t)]dr (Cl) 

^'->0 2 A a Jo 

where G is the energy release rate; A a is the crack extension at the crack tip; <t i2 (A a - r, 0) are 
the stresses prior to the crack extension; u, (r,± n) are the displacements due to extension. (Cl) 
can be rewritten by 

G = lim — f xl (A a - r, 0) [u{r,K) - u(r,-Tt)]dr (C2) 

Aa->0 2A a 

By Stroh formulation, the stress vector r 2 is written as 

r 2 (r.d) = *., = Re %B((8 m +l)z s ")B-'g m 

m= 1,2,3,- 

24 
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Along the crack plane, 0=0, 

T 2 (Aa -r,0) = + l)(Aa -r) 5 "'Re[g m ] 


m=I,2,3,~ 


= £(«„+lXA,» -!•)*■*„ 

= X O' + tXA a-ry~ yl g 2 ^ 


7 = 0 . 1 . 2 ,- 


Note that g m is pure imaginary for m = 2, 4, 6, — . For the displacement vector u, 

u = Re[ Y, A { zSn+ ') B ~' 


n= 1,2.3,- 


At the crack flanks, 9 = ±n , 


Thus 


/i_\ ±/ff 5 n -t-l 5„+l ±/(<5 n +1 )JT 

z = rcos(+;r) - re , z " = r " e ■ 


Am = u{r,n) - u(r,-n ) 

= Re{ £[A(r 5 " + V <5 " +1) * - 

/ i = l , 2 , 3 ,~ 

= Re{ Zr s ' + '2isin[(S„ +1 )n]AB^g n } 

11=1 . 2,3,- 

= 2 £ r 5 " +1 (-1) ( "- 1)/2 Re [IT 1 - iSU 1 ]g„ 

n=l,3,5, ••• 

= 2 ]£(-l)<^ )/2 r , s J1 / A < 

n = l , 3 , 5 .~- 

= 2 £(-l)' IT'*,,./ 


>1 


/+ 1/2 


/=0,1,2,~ 

-1 cr-l r :r-l 


where the identity - AB ~ 1 = SL 1 + /IT 1 has been used in deriving eq. (C4). 
Substituting (C3), (C4) into (C2) 


G = ^^^ 2U+ \ H - l) ' 8 ^' L " 8 -C /(Aa - r)J i^7 dr 


Introducing a nondimensional variable, x = r! Aa 


(C3) 


(C4) 


(C5) 


G = lisiSI 2 ( ^l )( - iy ^ J > ,(1 ~ xy ^& x dx 

In the integrand, the linear term of Aa corresponds to the first term orj = l = 0 


G = lim— !— [gjL~ l g,Aa-^ + o(Aa)] 
Aa^o 2 A a 2 

= -k T L~ ] k 


(C6) 


25 



For the monoclinic solid with plane of symmetry at xj, = 0, L' x under in-plane deformation can be 
expressed by 



b d 


d e 


(C7) 


/I, + fj. 2 = a + ib, fJ .^2 =c + id 
e = ad - be = Im [^fl 2 (/T, + p 2 )] 



b d' 
d e 



= sf , [ekf + 2 dk x k 2 + bk\ ] 

= ^,{lm[X// 2 (ft + p 2 )]k; + 2 Im(/ij// 2 )fc,fc 2 + Im(//, + // 2 )k 2 } 


(C8) 


This expression can be reduced to orthotropic solids without cross term of k\ and ki in eq. 
(C8). This special case has been formulated from Lekhnitskii formulation and the results are 
shown in eq. (26) of ref. (16) by G. C. Sih, P. C. Paris, and G. R. Irwin. 
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Appendix D 

Determining the Unknown Coefficients in the Crack Field Using the J-integral 

The path-independent integral J M from eq. (60) is 


r du a 

J M - J r (^,2 — r W,J ) ds 

(Dl) 

Let the actual and auxiliary crack fields 


n 

(D2) 

<P° =B{f a (z)}B-'h 

(D3) 


where 

f„(z) = Z S '+\S n =(n-2)/2, n= 1 , 2 , 3 ,-. 
f a (z) = z A ” + \ A m = —ml 2 , m= 1 , 3 , 4 ,-. 

From the identity 

Re(C) Re(D) = ~ Re[(C + C ) D] = ~ Re[C(D + D)] (D4) 

where C and D are complex matrices, we have 

< 7 j2 du° = du° Re[4>, , ] = d(u a ) T Re[d> , ] 

= cf[Re(AB _, O a ) 7 ']Re[0 1 ] (D5) 

= ^Re[(AB- 1 ^0") 7 (a> l +c.c.)] 

and 

d® a 

t"u t , = - Re( — )w ( , ds 

ds 

= -Re(dO a ) T Re[AB- ] O i ] (D6) 

= -jRe[(d©") 7 '(Air , ® 1 + c. c . )] 

where c.c. denotes the complex conjugate of the preceding term, i.e., 

F + c.c.= F + F 

Therefore 
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J„ =^Rejj r [(AB-'<i<t) , ) r l<l> l +c.c.) + (dQ> a ) T (AB~'&j + c.c.)]j 

=i Re |£J r ( h r B- T {dr)A T B(f:)B-g, + h T B- r (df)A r B{r;)B-'g , 

+ h T B r (df)B r A{f;)B-'g,+k r B 1 (dr)B T A(7;}B-'g,)} 

= 7 Re | 'L\S>‘ T B- T (dr)(A T B + B T A)(f;)B-'g„ 

+ h T B- T {df°)(A T B + B T A){J:}B- l g n )] 

Using the orthogonality relations in eq. (10) and (11), 

A t B + B t A = land A T B + B T A= 0, (D7) 

J M can be further rewritten as 

L " J -, (D8) 

= |Re J,h T B- T (\ r f n '(z)df a (zjjB- l g n 

Defining 

*_ - J r /:u«)<r u„) = <«. + ix a, + 1 

it is readily shown that 

(a) n = m or 8 n + A m + 1 = 0 , 

R m , = (S, + 1)(A„ + 1 )J ^ = -S,(S. + l)[ln(r 2 / r t ) + i2n ] ; 

J r 7 

(b) or <5 n + A m + 1 * 0 , 

= ( ^" ~~~ { [r?- m,n - - r["~ m)n ] cos(n -m)xf 2 + i{d n ~ m)n + r^ m)n ) sin (n -m)n! 2} 

(-1) ; ^ + 1)(2 ~ w) [r t«-mv 2 _ (—)/ 2 j for n-m = 2j 

n-m 

/ (- 1 y ( ^ +1)(2 " m) [r<"- m,/2 + r, ( ^ )/2 ] for «-m = 2y + l 

rc-m 

Using and the identity 

B t B =— 2i IT 1 

denoting 7*, as J v ,„„ /i as /i m , eq. (D8) becomes 
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J Mm ” 2 ^ 


X KXB- T B-'g„ 


= 2 Re 


J.KS-miL-'g, 


= Im "Y R„,„h T m L 'g. 

mn m o n 


= Im[fi,.*II.- | S .]+ I \m(R m „)hlL-'g,+ 2>(R„X£''ir.) 


The last two terms of the above equation are zero. Thus the term contributing to the J Mm is the 
term between g m and h m only, and 

J Mm = ^{Km h m T ^gn] 

= h m T L~ ] g m lm{-S m (8 m + l)[ln(r 2 !r x ) + Hk ] ) , n * 0 

= -2n5 m (5 m + Y)hJ L~' g m 


or 


J H ,—lnS,(S m *\)h, r L^g, 


(D9) 


As A n is arbitrary, we choose 


f e k , n = 1,3,5,- •• 
\ie k , n = 4,6,8, - 


Here, e^(k= 1,2) possesses dimension force / (length) 1-5 " . Therefore, for a given n, there will be 
two different values of J Mn denoted by J$ n and J { ^\ correspondingly. For the two choices of h n , 
eq. (D9) leads to 


where 


LJ 


M n 


gn =i 


2 k 8 n ( 8 n + 1 ) ’ 
i LJ Mn 

[ 2s:5,(5„ +1)’ 


n = 1 , 3 , 5 ,• 
n- 4 , 6 , 8 • 


t =u m y (2) f 

J Mn '- J Mn ' J Mn* 


(DIO) 


29 



Appendix E 

Determining the Unknown Coefficients in the Crack Field Using Betti’s Theorem 

Following arguments similar to those presented for the J M integral in Appendix D, one can 
get the expression for the integral I. The path-independent integral from Betti’s theorem is 

l = \ r [{u a ) T t-(t a ) T u]ds (El) 

Using the complex potential functions, O and <E> a and identity eq. (D4), 

(u a ) tds = - Re( AB O a ) T ] Re^ZO] 

1 (E2) 

= --R e[(AB- ] d® a +c.c.) T d®] 

(t a ) T u ds = -Re(d<t> a ) T Re[AB-'®] 

= -^Re{[J(O n +c.c.) t }AB~'<1>} (E3) 

= -^Re{J[(O a + c.c.) r A2? _I <J>] - (<t> a +c.c.) T AB~'d<&} 

j (f 2 , X) 

j r (t a ) J uds = --Re{[(®“ +c.c.) t AB- 1 ®] -J (<J>° +c.c.) r AB"VO} (E4) 

According to eq. (19), the traction free conditions on the crack faces for the auxiliary field may 
be written as 

Re[O a (z)] = 0 or <J>° +c.c.= 0 at <9 = ±tt . 

Therefore the first term on the right hand side of (B4) vanishes and 

J r (C) T M J5 = |Re[J r (<D° +c.c.) r AB‘'#] (E5) 

Substituting (E2) and (E5) into (El) and using the expressions for 0 and O a (eq. (D2) and (D3)) 
yield 

/ =-iRejj r (AB'<J>“ +c.c .) 7 

= -T Re |Xj, h T B- T (f)A T B(df,)B-'g, +h r B- T {f°)A T B{df,)B-'g, 

+ h T B- T {f)B T A{df,)B-'g,Ah T B- T (f)B T A{df,)B-'g „] 

/ = -(-Re XL li T B- T (r)(A'B + B r A,tdf.)B-'g, : + 

^ L n 

[h r B- T (ryA T B + B T A)(df,)B-'g,]] 
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Using the orthogonality relations, eq. (D7), 

' = 4«4lJ r h T B- T {r){df,)B-'g, 

^ L n 

= -yRe Xft r B' r {j r /‘(z)#,(z))B''s. (E6) 

where 

f a (z) = z K+l , A m = -ml 2 , m = 1,3,4,-. 

Defining 

&„ - J r r (z. w, (z„ )=<«,+ 1 ) J r 

it is readily shown that 

(a) « = m - 2 or 5 n + A m + 2 = 0 

Qm.m- 2 = (d m -2 + ^)J f = (^m-2 + l)(^ nr 2 ^ K i + , 

(b) n^m -2 or 8 n + A m + 2 * 0 

Q„„ = -2 - r/"^ 2 ’ 72 ) cos[(n - m)K / 2] + i(r^~ m+2)n - +r 1 '”"" +2),2 )sin[(n-mk/2]} , 

(_ d 2 J^L±lL ( r («- m ^)/2 _ (»-«* 2)« ) for n-7n = 2; 

n-m+2 ‘ 

2 i (- 1) 2+1 — + - (r, (fl ~ m+2)/2 + r 1 (n ~ m+2)/2 ) for n - m = 2j + 1 

n-m+2 ~ 

Using 2™ and expression for L’ 1 , writing 7 as /„ , and h as 7t m , the 7 -integral from eq. (E6), for a 
given A m , becomes 

L = -\mY,Q.XB' T B-'g,} 

** n 

= -|Re(Xe,.(- 2 o*:fs.i = -lincXa. KB~'g,i 

=-ta>(e„,.. 2 *iL-' f „. 2 )-im Xi&„ />:£''*„ i-i”> 

n n 

y , /i<«-2 m=m+2j + l 

Since the last two terms of the above equation are equal to zero, it can be clearly seen that the 
term contributing to the 7 is the cross terms between g m .2 and h m only. 

Thus 

7 m = ~1n{Q m , m .2K T L~'g m -7 ] = ~K T L~'g m - 2 Im{(<5„,_ 2 + DBn (r 2 / r, ) + Hk}} 

= -2n( 8 m _ 2 + \)h m T L~ x g m _ 2 
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° r '..i=-2’t(5 n + l)h n ,/L- , g, 

From eq. (E7), following a similar procedure as before, g n can be expressed by 

LI „+ 2 , t c 

-zm — rr» « = i,3,5, 

J 2 x( 8 „+l) 

Sn } iLI n+1 

r- 5 ^ — , n = 2,4, 6,-- 

2k (5 n + 1 ) 

where 


I =r/ (n / (2) l r 

* n+2 ’ * n+ 2-1 


/'*2 is the value of 7 n + 2 when 


e* , n = 1 , 3 , 5 ■ 

ie. , « = 2, 4, 6,-- 


Here, possesses dimension force x (length) 5 " . 
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